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Recent progress in the calculation of the two-loop on-shell mass counterterms within the electroweak Standard 
Model (SM) for the massive particles are discussed. We are in progress of developing a package for full two- 
loop SM calculations of 2 — > 2 fermion processes, with emphasis on the analytical approach where feasible. The 
complete two-loop on-shell renormalization is implemented. Substantial progress has been made in calculating 
the master integrals. We are able to compute in an efficient and stable manner up to a few thousands of diagrams 
of very complex mass structure. 



1. Introduction 

Even for the simplest physical processes like 
e+e^ ff complete two-loop electroweak SM 
calculations still are not available, maninly be- 
cause of the enourmous complexity of such cal- 
culations. An exception is the /i-decay rate 0, 
which is simpler due to the fact that it is a static 
quantity. Besides the large number of diagrams 
encountered, the difficulties start at the level of 
individual Feymnan integrals and increase sub- 
stantially when going from propagators to ver- 
tices or box contributions. One key problem at 
the beginning is the renormalization and the cal- 
culation of the necessary counterterms. In the 
QED like on-shell renormalization scheme the ba- 
sic input parameters for electroweak higher or- 
der calculations are the fine structure constant 
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a = and all the physical particle masses. In 
this scheme the whole renormalization program 
only requires the calculation of selfenergy dia- 
grams and tadpoles. While the calculation of 
the counterterm for the fine structure constant 
is relatively easy (at zero momentum) [2j, the 
on-shell mass counterterms are much more in- 
volved. In this note we therefore focus on as- 
pects of calculating the latter, which yields the 
relation between bare, MS and on-shell (pole) 
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masses (two-loop renormalization constant in on- 
shell scheme). Many details on this program may 
be found in the original publications j3»4j and will 
not be repeated here. We rather concentrate on 
some controversal points concerning the defini- 
tion of MS masses in electroweak theory (Sec. 2). 
We also give some more specific information con- 
cerning our numerical approach to the calculation 
of huge sets of diagrams (Sec. 3). 

2. MS-masses of particles in SM and renor- 
malization group equations 

The two-loop calculation of pole masses of 
the gauge bosons in the SM has been discussed 
in |3I4| . In terms of the transversal self-energy 
function , vn? , • • •), by expansion about = 
— m^, will get the two loop solution 

sp = m2_n(i)_n(2)-n(i)n(i)', (i) 

for the location of the pole sp. 11^^^ is the bare 
or MS-renormalized L-loop contribution to 11, the 
prime denotes the derivative with respect to p^. 
One of the remarkable properties of 0J is that the 
complex pole is represented by the self-energy and 
its derivative at momentum equal to the bare or 
MS mass which, by definition, are real parame- 
ters. The standard parametrization of the pole is 
sp,a = — iMaTa, whcrc Ma is the pole mass 
and Fa is the width of particle a. 

For calculations of electroweak corrections in 
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the SM two renormalization schemes are com- 
monly accepted: the on-sheU and the MS scheme. 
It is well understood, that in the on-shell scheme 
all momentum independent diagrams, in partic- 
ular the tadpoles, can be omitted (the on-shell 
scheme is a particular case of a momentum sub- 
traction scheme: finite parts are fixed by subtrac- 
tion of the propagators at = —sp). The set 
of diagrams contributing to the MS mass are not 
unambiguously defined. Let us first express the 
pole 1^ in terms of the bare amplitude in a man- 
ifestly gauge invariant manner. This requires to 
include the Higgs tadpole contribution 5 . Only 
this complete gauge invariant bare amplitude 
should be utilized as a starting point to set up 
MS renormalization. At the two-loop level MS 
renormalization can be written as 



Sp = rriQ 
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where the sum runs over all species of particles, 
gj = a, gs, (<55j,o)^^^ and {Sm'j QY^'> are the one- 
loop counterterms for the charges and physical 
masses in the MS -scheme and after differentia- 
tion we put all parameters equal to their on-shell 
values. The derivatives in Eq. correspond to 
the subtraction of sub-divergencies. The genuine 
two-loop mass counterterm comes from the shift 
of the rriQ term. The relation between bare- and 
MS -masses has the form 



k=l 



(3) 



To renormalize the pole mass at the two-loop 
level requires to calculate the one-loop renor- 
malization constants for all physical parameters 
(charge and masses), and the two- loop renormal- 
ization constant only for the mass itself. Not 
needed are the wave-function renormalization or 
ghost (unphysical) sector renormalizations. Af- 
ter UV-renormalization the pole is represented 
in terms of finite amplitudes. Now, expression 
^ connects the pole sp with the MS parame- 
ters: masses and charges. This expression can be 



inverted and solved iteratively. The solution to 
two-loop reads 

m^^Mj+R.{n«}_+Re{n?>+n>'>n<'>'}_ 



(A«)")| + 2(Am3)'"^ 
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where the sum runs over all species of particles 
J = Z, W, H, t, (Am^OW = Re{n,}^ , and 
the transition from the MS to the on-shell scheme 
for the electric charge [5] is also included. The 
mass on the l.h.s. of this expression we call the 
MS-mass of particle. It should be noted, that in 
this definition the tadpole contribution does not 
cancel, so that higher powers of the Higgs and 
the top-quark mass show up at higher orders. In 
particular, at two-loops, the purely bosonic di- 
agrams generate m|^/my terms and the third 
fermion family gives rise to the appearance of 
mf/{m'jjmy) power corrections. For the MS- 
masses, defined in this way, the following prop- 
erties are valid: 

1. The UV counter-terms satisfy relations con- 
necting the higher order poles with the lower or- 
der ones: 



la' 



d 



dm?: 



1 Vo — 

2 2^33 Qg^Z a , 



(5) 



where we adopt the following definitions for the 
RG functions: for all dimensionless coupling con- 
stants, like g, g' , gs, e. A, yt, the /3-function is given 
by fJ^-^i^g = Pg and for all mass parameters (a 
mass or the Higgs v.e.v. v) the anomalous dimen- 



RG-invariant: 

^ sp = 0, we are able to calculate the anoma- 



sion is given by /i In m - 
2. Using the fact that sp is 

lous dimension of the masses from our finite 
results (0} or from the UV counterterms © 
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3. All tree level relations between masses of any 
particles and parameters of the unbroken La- 
grangian are RG invariant. This means, in par- 
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ticular, that the RG equation for the vacuum ex- 
pectation value V is given by 7^2 = 7^2 — P\/X, 
where and A are the parameters of the sym- 
metric scalar potential. This fact allow to get 
anomalous dimension of the masses via the rela- 
tions [3 ,j 

7w = 7m2 - — -I- 2— , 
A 9 

/3a Pg^ f3g'\ 

A \ a 9 ) 

/^A . f^yt ia\ 
7t = 7m2 - H , IH = 7m2 , (6) 

A yt 

where sw{cw) are the sin (cos) of the weak 
mixing angle and the 2-loop RG functions 
Pg, Pg' 1 Px,lm^T Pyt ^^rc Calculated in the unbro- 
ken phase 0. 

The RG invariance of the pole positions sp allow 
us to factorize explicitely the RG logarithms 

where Lb = In^^/m^, c["^:^^ = cj™;"^ and 

^^i^^-^=l''=Y.a^lf^9i9^l'^ 

j 1,3 

d 

^^•^^ ^J-^ + i (J-)Jb) ib) ia)\ .J. . ^ 

2 ^2 y^t ^^2 3 13 ^^2 i y 

^cSS-^P3ifk3+^'r'if 

,^ f (*-■) 2 ^ (a) , (fc) 2 (a)A 

where X^"''^ and their derivatives can be extracted 
from Appendixes of 5 3 4. 

Crucial point of our definition of the MS -mass 
Q is the gauge invariant construction ^ for the 
pole in terms of the unrenormalizcd, bare dia- 
grams. It can be done only after inclusion of the 



Higgs tadpole contribution. Another important 
ingredient are the Ward identities. 

A. The inclusion of the tadpoles is necessary to 
ensure, that the physical Higgs field has zero vac- 
uum expectation value in each order of the loop 
expansion. 

B. It is well know, that in order to preserve the 
Ward identities for the longitudinal part of the 
gauge boson propagator it is necessary to add the 
tadpole contribution, which is equal to the prop- 
agator of the would-be- Goldstonc bosons at zero 
momentum transfer. In particular, at the two- 
loop level, the photon would aquirc a mass if the 
tadpole contribution would be omitted. 

Our RG equations for the v.e.v. v and the 
particle masses m are different from the ones ob- 
tained in the effective potential approach A 
comparison of predictions based on these two ap- 
proaches have been recently performed in IIO'. 

These structural considerations were important 
to check our calculations of the various counter- 
terms. 

For the calculation of the 0{aas) and the 
O(a^) corrections to the top-quark propagator we 
refer to @] and respectively. 

3. Numerical results 

According to we need to calculate 

propagator-type diagrams up to two loops on- 
shell. To keep control of gauge invariance we 
adopt the gauge with three different gauge 
parameters ^w, £,z and For our calculation 
all diagrams have been generated with the help 
of QGRAF U^. The C-program DIANA 
then was used together with the set of Feynman 
rules extracted from the package TLAMM 
to produce the FORM input which is suitable 
for the package ONSHELL2 15 and/or for 
another package based on Tarasov's recurrence 
relations ^Hl- The set of master- integrals, in the 
limit of massless lepton and light quarks (see de- 
tails in 171) includes diagrams with three different 
massive scales. In most cases exact analytic re- 
sults in terms of known functions are not avail- 
able. Thus, instead of working with the exact 
formulae (which only can be evaluated numeri- 
cally, at present) we resort to some approxima- 
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tions, namely, we perform appropriate series ex- 
pansions in (small) mass ratios. For diagrams 
with several different masses it is possible that 
several small parameters are available. In this 
case we apply different asymptotic expansions 
(see ^Tj) one after the other. Specifically, we 
expand in the gauge parameters about = 1, 
in siia^ 9w and, for diagrams with Higgs or/and 
top-quark lines, in rriy/m'jj or/and my /ml. Nu- 
merical results are obtained using the packages 
ON-SHELL2 15 and TLAMM [Tl|. Since 
the quality of the convergence of a series is not 
known a priori we have to calculate several coef- 
ficients of each expansion (six in sin^ Ow and five 
in mass ratios, my / {mfj , m'^) ) to keep control 
on the convergence. For the one-loop diagrams 
and their derivatives we used the exact analyti- 
cal results, as given in |18i;-{| . The expansion of 
diagrams with a top-quark and/or a Higgs bo- 
son leads to two-loop bubble diagrams with three 
massive lines (with two of the masses equal). For 
these master integrals we utilized a special form 
of representation 18 3]- The diagrams with mass- 
less fermion lines also demand special considera- 
tion. These diagrams develop threshold singular- 
ities which behave like powers of Insin^ 9w- To 
control these terms we had to use the exact ana- 
lytical results, which have been worked out in P) 
using a technique developed in |HI18| . We found 
that after collecting the contributions from all di- 
agrams the threshold singularities canceled. This 
is a manifestation of the infrared stability of the 
pole mass of the gauge bosons. The series ex- 
pansion in sin^ 9w converges very well, and can 
be restricted to the first two coefficients. The 
expansion in the remaining mass ratios require 
three coefficients in order to get sufficient preci- 
sion for light Higgs mass values. The numerics 
has been performed in MAPLE. To get control of 
the numerical stability, we run the MAPLE pro- 
gram with an accuracy of 100 decimals (a pos- 
teriori, as an experimental fact, we find that the 
minimal accuracy is 40 decimals). 
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